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ABSTRACT 

A method is presented for improving the numerical prediction of bypass transition heat transfer on a flat plate in a 
high-disturbance environment with zero or favorable pressure gradient. The method utilizes low Reynolds number 
k-e turbulence models in combination with the characteristic parameters of the transition region. The parameters 
representing the characteristics of the transition region used are the intermittency, transition length and turbulent 
spot properties. An analysis is made of the transition length in terms of turbulent spot variables. The nondimensional 
spot formation rate, required for the prediction of the transition length, is shown by the analysis to be a function of 
the spot spreading angle, the dimensionless spot velocity ratio and the dimensionless spot area ratio. The 
intermittency form of the k-e equations were derived from conditionally averaged equations which have been shown 
to be an improvement over global-time-averaged equations for the numerical calculation of the transition region. 

The numerical predictions are in general good agreement with the experimental data and indicate the potential use of 
the method in accelerating flows. Turbulence models of the k-e type are known to underpredict the transition length. 
The present work demonstrates how incorporating transition region characteristics improves the ability of two- 
equation turbulence models to simulate bypass transition for flat plates with potential application to turbine vanes 
and blades. 


INTRODUCTION 

The present effort is part of a program by the National Aeronautics and Space Administration (NASA) for the 
investigation of the heat transfer in the transition region of turbine vanes and blades. 

The transition from laminar to turbulent boundary layer flow affects great increases in the local wall shear stress and 
heat transfer. This effect of transition is especially critical for airfoil surfaces such as turbine blades/vanes where 50 
percent or more of the surface can be in transition. The ability to predict the onset and extent of transition is impor- 
tant to the determination of turbine blade heat transfer that critically affects durability and engine thermal perfor- 
mance. The accurate prediction of gas side shear stress and heat transfer on turbine blades will become more critical 
as the desired operating temperatures of advanced turbine engines increase. The present methods for predicting 
transition shear stress and heat transfer on turbine blades are based on incomplete knowledge. There is incomplete 
knowledge about the transition process in an engine environment where disturbance levels are initially large, 
although recent surveys (Mayle, 1991; and Volino and Simon, 1991) indicate much progress in this area. In these 
large disturbance environments traditional linear mechanisms are bypassed and finite nonlinear effects must be 
considered. 

Full modeling of this transition phenomena for the purpose of predicting wall shear stress and heat transfer on 
turbine blades must consider the effects of free-stream turbulence, pressure gradient, streamwise curvature, surface 
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roughness, wall/free-stream temperature ratio and flow disturbances (e.g., wakes). The present work, a first step, 
concentrates on the effects of free-stream turbulence and pressure gradient on a flat plate. The objective is to 
demonstrate a modeling approach to the transition region that will improve the numerical prediction of that region. 

Properly used. Low Reynolds Number k-£ turbulence models for calculating the transition heat transfer in environ- 
ments of high disturbance appear to simulate the onset of transition governed by the transport of turbulence energy 
from the free stream into the boundary layer. However, while two-equation turbulence models predict an acceptable 
value for transition onset, the resulting process is generally too abrupt, resulting in an underprediction of the 
transition length (e.g., Rodi and Scheuerer, 1985; and Simon, 1993). The experimental evidence (Volino and Simon, 
1993) suggests that the physics of the transition region is characterized by an incomplete development of the cascade 
of energy from large to small scale turbulence and the k-£ model does not completely account for this. Three 
modifications of the k-e model, as a potential means of improving the simulation of transition physics, are reviewed 
by Simon (1993). The three methods are (1) the use of a multi-time-scale k-e model (Crawford, 1992), (2) modifica- 
tion of the turbulent production term in the differential equation for the turbulent kinetic energy (Schmidt and 
Patankar, 1988), and (3) the use of intermittency in the k-e equations (Simon and Stephens, 1991). 

Simon and Stephens combined the conditionally averaged momentum, energy, and turbulence equations for the 
turbulent spot and nonturbulent portion of the intermittent flow of the transition region to obtain global values of 
velocity and temperature. The resulting equations contain the additional variable of intermittency which requires 
knowledge of the transition length. A transition length estimate was made by Simon and Stephens (1991) using the 
method of Narasimha (1985). It expresses the transition length in terms of a nondimensional spot formation rate (N) 
which is assumed to be constant for free-stream turbulence levels greater than about 1 percent. Simon and Stephens 
followed the approach of Narasimha and determined the value of N based on experimental data. This permitted a 
calculation of transition length in terms of the momentum Reynolds number for transition onset. The transition 
length may also be determined using an empirical approach (e.g., Dhawan and Narasimha, 1958; and Narasimha, 
1978) or by the use of turbulent spot analyses (Walker and Gostelow, 1990; Walker, 1989; McCormick, 1968). 

The present work follows the line of reasoning given in the latter analytical works to develop a model for the 
nondimensional spot formation rate (N) for a prediction of transition length and intermittencies required for the 
equations of Simon and Stephens. The equation developed for predicting of the transition length is expressed in 
terms of turbulent spot characteristics of spreading angle, velocity, and physical dimensions. A check of the analyti- 
cal formulation is made using the recent turbulent spot data of Clark, LaGraff, Magari, and Jones (1992). 

Having a method to predict transition length and intermittency permits the prediction of transition heat transfer. In 
this work numerical predictions are made for transition heat transfer, using the modified Jones-Launder turbulence 
model of Simon and Stephens (1991) with the TEXSTAN boundary layer computer code (Crawford, 1985). The 
numerical predictions are compared with the flat plate, zero, and favorable pressure gradient data of Blair and Werle 
(1980,1981). 
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constants appearing in the k-e turbulence model 
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transition length 
Mach number 

nondimensional spot formation rate 
spot formation rate 
static pressure 
turbulent Prandtl number 
wall heat flux 

momentum thickness Reynolds number 

Stanton number 

temperature 

fluctuating temperature 

turbulence intensity at free-stream 

turbulent heat flux 

free-stream velocity 

turbulent shear stress 

mean velocity in x,y,z directions 

fluctuating velocities in x,y,z directions 

location of transition onset (spot formation begins) 

spreading angle of spot 

velocity of center of spot/free-stream velocity 

intermittency 

dissipation rate 

boundary layer momentum thickness 
area of spot/square of half width 
thermal conductivity/specific heat 
molecular viscosity 
eddy or turbulent viscosity 
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V 


kinematic viscosity 


p fluid density 

a dimensionless spot propagation parameter 

a k , o e empirical constants in turbulence model 

Subscripts 

e free-stream value 

tr transition onset 

E transition end 


Transition Region Equations 

The beginning of the transition region is defined by the first appearance of turbulence spots. The transition region is 
characterized by the intermittent appearance of these turbulent spots, which grow as they move downstream until 
they finally merge to form the turbulent boundary layer. The transition region, therefore, may be described as having 
turbulent spots with laminar-like fluid surrounding the spots. Vancoille and Dick (1988) state that conventional 
turbulence models based on global time averages cannot give a good description of such intermittent flow, resulting 
in poor agreements for the turbulence intensity profiles for the boundary layer. They reported that using conditional 
averaging techniques for the turbulent spots and the laminar like fluid surrounding the spots resulted in a good 
prediction of the experimental values of the velocity profiles and the turbulence intensity profiles in the transition 
region. The work of Simon and Stephens (1991) extended the work of Vancoille and Dick (1988) to predicting heat 
transfer in the transition region. The conditioned momentum, energy and turbulence equations for the turbulent spots 
and nonturbulent portion (“laminar”) of the intermittent flow in the transition region are combined and simplified to 
obtain global values of velocity and temperature. The present work refines the previous work of Simon and Stephens 
(1991) by removing the need to use an empirical nondimensional spot formation rate (N) for prediction of transition 
length and by extending the prediction of N to include pressure gradient. A brief review of the key equations follow. 

As shown by Simon and Stephens, when the modeled equations for the turbulent and laminar zones are combined 
and the assumption made of negligible Reynolds stresses in the “laminar” zone, the following momentum and 
energy equations are obtained: 
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Equations (1) to (6) are the simplified equations for the transition region derived by Simon and Stephens (1991). 
Constants for the turbulence model are given in Table I. These equations require the specification of intermittency. 
Equation (1) shows the Reynolds stresses multiplied by the intermittency factor. This is the form utilized by 
McDonald and Fish (1973) and others. 


TABLE 1— CONSTANTS AND FUNCTIONS FOR JONES-LAUNDER 
TURBULENCE MODEL 


n 

Ci 

C 2 


(a) k-c Turbulence model constants 


0.09 

1.45 

2.0 

1.0 

1.3 


(b) Functions 


f u exp [-2. 5/(1 + R-^50)] 
f 2 1.0 -0.3 exp -R^ 

R t pul / jx — pk 2 /ep 


Intermittency 

Specification of intermittency requires knowledge of the transition path in terms of transition start and of the 
transition length. Narasimha (1957), utilizing a concentrated breakdown hypothesis, derived the following transition 
path equation from the turbulent spot theory of Emmons (1951): 
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Equation (7) has been compared to experimental data and has been found to be valid. Figure 1 , from Volino and 
Simon (1991), compares the experimental data from Kim (1991), Kuan and Wang (1990), Sohn, O Brien, and 
Reshotko (1989), and Blair and Anderson (1987) with Eq. (7). Figure 1 is an example of the good agreement that has 
been found with experimental data and Eq. (7). 

Using the approach of Narasimha (1985) the transition length may be expressed in terms of the transition Reynolds 
number and a nondimensional spot formation rate as follows: 
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Figure 1. — Intermittency (Volino and Simon, 1991). 


Where N, the nondimensional spot formation rate is defined as 

N = ncr0^/v (9) 

In what follows an analysis is made of spot inception and growth. Spot inception is described in terms of the onset 
position and frequency of disturbance. The frequency of disturbance is taken to be the Tollmien-Schlichting wave 
frequency having the maximum amplification. This permits an expression for the dimensionless spot formation rate 
(N) in terms of the spot spreading angle, spot velocity, and spot physical dimensions. The transition length may then 
be calculated using Eq. (8). 


Spot Formation Rate (N) Analysis 

Chen and Thyson (1971) used the observations obtained from stability experiments with controlled disturbances, 
that prior to laminar flow breakdown the wave motions become three-dimensional with regularly spaced peaks and 
valleys. Breakdown or the appearance of turbulent spots is assumed to occur at the peaks. The turbulent spots grow 
and spread and eventually meet to form a fully developed (albeit immature) turbulent boundary layer. This process is 
illustrated in Fig. 2. The transition length, defined as the length from inception to the first meeting of adjacent spots, 
was shown by McCormick (1968) to be related to the distance between turbulent spots at onset (see Fig. 2) as 
follows: 

d - 21^ tana (10) 

From the following intermittency-path equation derived by Narasimha (1957) 
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Substituting Eq. (12) in Eq. (10) we obtain 


U 

d = 2| 4.61 — - 
n<r 


1/2 


tana 


(13) 


Based on the assumption that breakdown to a turbulent spot occurs at the peak for each cycle of wave motion, Chen 
and Thyson (1971) state that the spot formation rate (n) can be related to the frequency of disturbance (f) and the 
spanwise wavelength. This is stated here in terms of the distance between onset of spots (d) as follows: 


n = f Id 


(14) 


So that by combining Eqs. (13) and (14) we have 


4tan 2 a4.61U e 

The dimensionless spot formation rate as defined by Narasimha (1985) is given by Eq. (9), which can be reformu- 
lated and expressed as 


N = 


nov 2 




Substituting Eq. (15) in Eq. (16) we obtain 
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(17) 


According to Eq. (17) the dimensionless spot formation rate is a function of disturbance frequency, dimensionless 
spot propagation parameter, the turbulent spot spreading angle, the momentum Reynolds number for transition onset 
and the free-stream velocity. For the purpose of estimating N, the frequency (f) of spot formation is determined by 
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using the method of Walker (1989). Walker applied the model of McCormick (1968) who assumed that turbulent 
spots appeared at the frequency of the basic Tollmien-Schlichting instability wave. He assumed that the dominant 
disturbance frequency at breakdown to be the TS wave frequency having the maximum amplification ratio accord- 
ing to the small-disturbance stability calculations, for a Falkner-Skan profile, obtained by Obremski et al. (1969). 
The equation for the locus of maximum amplification rates over the whole range of the Falkner-Skan pressure 
gradient parameter (including zero pressure gradient) as given by Walker is as follows: 


0.7676Ug 
2tcvR el 12 


(18) 


Walker presents experimental evidence, for a free-stream turbulence of 0.3 percent, to justify employment of Eq. 
(18). Roach and Brierley (1992) also demonstrate the ability of Eq. (18) to predict maximum T-S frequencies for 
transition occurring at 1 percent free-stream turbulence and zero pressure gradient. Substitution of Eq. (18) in Eq. 
(17) results in the following 


N = 8.1x10“ 4 [— ^ — 

Vtana 

From Emmons (1951), the following expression for the dimensionless spot propagation parameter is obtained 
A tan 2 a 


(19) 


( 20 ) 


Where A is a dimensionless parameter related to the physical dimensions of the spot and 6 is a dimensionless spot 
velocity. 

Substituting Eq. (20) in Eq. (19) the following equation for N is obtained 


N = 


f A tan a N 

l - J~) 


( 21 ) 


Equation (21) for N is only a function of turbulent spot characteristics and appears to justify the use of a constant 
value for N as suggested by Narasimha (1985). Narasimha (1985) found that for freestream turbulence levels greater 
than 0.1 percent the value of N has the approximate constant value of 0.7x1 0" 3 . However, using an 1 1 ° spot spread- 
ing angle as determined by Schubauer and Klebanoff (1955), a value of A = 2 as estimated by Emmons (1951) and a 
value of 6 = 0.65 based on the work of Clark, Jones, Ashworth, and LaGraff (1991) results in a value of N = 
0.29xl0' 3 from Eq. (21). According to Eq. (20) the value of a is 0.12. On the other hand, Narasimha (1985) reports a 
G value of 0.25 which, using Eq. (19), results in N = 1.3xlO' 3 . Thus we can see that N is quite sensitive and there is 
a fair spread in its value. However, the calculation sensitivity is lower for the transition length since it is a function 
of the square root of N (Eq. (8)). Using the experimental data sets of Schbauer and Skramstad (1948), Abu-Ghannam 
and Shaw (1980), Sohn and Reshotko (1991), Suder, O’Brien, and Reshotko (1988), Kim (1991), and Eq. (8), values 
of N were calculated and are given in Fig. 3. In Fig. 3, the values of N determined from the experiments are com- 
pared to the above calculated values of N using Eqs. (19) and (21). As shown in Fig. 3 the calculated N values cover 
the range of the experimental data, with a preference for the lower N value of 0.29x1 0' 3 . In the above calculations 
for N at zero pressure gradient, a spreading angle of 1 1 ° was assumed for the low Mach number conditions of the 
experiments. Clark, Jones, and LaGraff (1994) report that the spot spreading angle decreases monotonically with an 
increase in the Mach number. Their experimental data is shown in Fig. 4. Extrapolation of their data to the low Mach 
number conditions of the experiments of Fig. 3 indicates that the use of an 11° spreading angle is justified. 
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The experimental turbulent spot data reported by Clark, LaGraff, Magari, and Jones (1992) can be used as an 
example of the utility of Eqs. (8) and (21) to calculate the transition length. The experiments presented in this paper 
were performed at Oxford University in a 6-in. Isentropic Light Piston Tunnel. Measurements were made of 
naturally occurring turbulent spots with the use of surface thin-films heat transfer gauges. Assuming A = 2 and B = 
0.65, calculations were made for the transition length and compared to experimental values obtained from 
intermittency plots presented in the above reference. The results are shown in Table II for a zero pressure gradient 
and a mild adverse pressure gradient. The comparisons between experiment and calculation in Table II are good and 
indicate that the method presented in this paper shows promise. 

Utilizing the value of N = 0.29x 1 O' 3 obtained above in Eq. (8) results in the following equation for the calculation of 
the transition length for a zero pressure gradient: 

Re. =124Rep /2 for y = 0 -> 0.99 ( 22 > 

L lr 0 tr 

Wyganaski (1981) reports a turbulent spot spreading angle of 5° for a favorable pressure gradient. The transition 
length equation, based on a 5° spreading angle, for a favorable pressure gradient with a k = 0.75x1 O' 6 , assuming the 


□ Schbauer and Skramstad (1948) 

O Abu-Ghannam and Shaw (1 980) 

O Sohn and Reshotko (1991) 

D Suder, O'Brien and Reshotko (1988) 
A Kim (1991) 



Figure 3. — Variation of nondimensional spot formation rate with 
free-stream turbulence (Narasimha (1 985) method). 



Figure 4. — Variation of spot spreading angle (a) with free-stream 
Mach number (MJ for zero pressure gradient. 
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TABLE EL— COMPARISON OF CALCULATED TRANSITION LENGTH WITH THE 
EXPERIMENTAL DATA OF CLARK, LaGRAFF, MAGAR., AND JONES (1992) 


Pressure 

gradient 

a, 

Experimental 

N, 

Eq. (21) 

M ft ) 




Calculated Experimental 

Eq- (8) 

Zero, 

M e = 0.55 

7.6° 

0.137 x 10' 3 

0.28 

0.20 

Mild adverse 

15° 

0.55 xlO 3 

0.14 

0.16 


same A value used in the zero pressure gradient case and a 6 = 0.8 reported by Wyganaski, is as follows: 

Re L = 344 Re^ 72 for k = 0.75xl0“ 6 (23) 

Equations (22) and (23) may be graphically represented in the manner suggested by Dhawan and Narasimha (1958). 
This is done in Fig. 5 for a number of experimental data sets given in a survey report by Volino and Simon (1991). 
Figure 5, for the cases of flow acceleration, uses the average velocity between transition onset and end for calculat- 
ing the Reynolds number for the transition length (ordinate of Fig. 5). In general, there is good agreement of Eqs. 

(22) and (23) with the experimental data. From the pressure gradient conditions reported by Wyganaski (1981), a 
value of k = 0.63x1 O' 6 was calculated at onset for the Blair and Anderson (1987) experiment (k = 0.75x1 O' 6 ). This 
helps to justify the use of a 5° spreading angle for obtaining Eq. (23). Figure 5 demonstrates that with an increase in 
flow acceleration there is an apparent increase in transition length and that this increase in transition length is 
consistent with the characteristics of turbulent spots under accelerating conditions. The above is not considered to be 
a definitive test for the use of Eq. (21), but points to the need for a further study of turbulent spots subject to 
conditions of pressure gradient. 
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Figure 5. — Transition length correlation. 
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Heat Transfer Calculations 


A calculation of the transition region was made utilizing Eqs. (1) to (7) and Eqs. (22) and (23) for estimating the 
transition length. These equations were numerically solved by using the TEXSTAN boundary layer computer code 
(Crawford, 1985). TEXSTAN is based on the STAN5 boundary layer program developed by Crawford and Kays 
(1976). The finite difference scheme of TEXSTAN is based on the numerical algorithm by Patankar and Spalding 
(1970). TEXSTAN solves the steady two-dimensional parabolic differential equations that govern boundary layer 
flow. This program sequentially solves the momentum equation and any number of transport equations such as 
stagnation enthalpy, turbulent kinetic energy (TKE), turbulent dissipation rate (TDR), and mass concentration 
governing equations. A calculation for the intermittency was initiated at transition onset. Transition onset was 
determined by the numerical method employed by Simon and Stephens (1991). This was assumed to occur when 
the boundary layer calculations prior to the use of intermittency showed a rapid increase in turbulence kinetic energy 
with a resulting increase in skin friction. Predicted onset values are compared, in Table III, to the experimental 
correlations of Abu-Ghannam and Shaw (1980) and Mayle (1991 ). Mayle’s correlation does not account for pressure 
gradient; however, the effect of pressure gradient is small for a k of 0.75x10 In general, there is a good compari- 
son between the predicted and experimental transition onset values. 


TABLE III —COMPARISON OF PREDICTED TRANSITION MOMENTUM THICKNESS REYNOLDS 
NUMBERS WITH EXPERIMENTAL CORRELATIONS 


T , 

U’ 

percent 

Pressure 

gradient 


Re fi 

U tr 


Predicted 

Experimental correlations 

Method of 
Simon and 
Stephens, 1991 

Abu-Ghannam 
and Shaw, 
1980 

Mayle, 

1991 

1.4 

Zero 

356 

410 

324 

2.8 

Zero 

236 

210 

224 

2.2 

k = 0.75xl0' 6 

244 

282 

244 


The following initial and boundary conditions were used in the numerical calculations. 


Initial and Boundary Conditions 

Wall boundary condition. Along the wall the no-slip boundary condition was applied. 


y = 0 : u = v = 0 


For the energy equation a wall heat flux was prescribed. 

q w=- K |r ( y =0) 


(24) 


(25) 


The turbulence kinetic energy and dissipation rate are set to zero at the wall, 
y = 0 : k = £ = 0 

The zero boundary condition for the dissipation was made possible by Jones and Launder (1973) by including 
additional terms in the turbulence energy (Eq. (5)). 
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Freestream boundary conditions. The freestream conditions may be expressed as a set of coupled ordinary differen- 
tial equations as follows: 


U 


e 



(27) 



(28) 


Where initial values of the free-stream turbulence kinetic energy are determined from the turbulence intensity as 
follows: 


k e='- 5 ( Tu e U e) 2 (29) 

Jl/3p + ^+^) 

where Tu e = — (30) 


And the initial dissipation rate determined from the following: 

£ e= k e /2/L e (3 0 

The dissipation length scale (Le) can be calculated from the longitudinal integral length scale using the method of 
Hancock (1980). 

Initial values. To initiate calculations, a Blasius velocity profile and a flat temperature profile (unheated start) were 
utilized. For the turbulence transport equations the following initial profiles suggested by Rodi and Scheuerer (1985) 
were used: 


k = k e( u/U e) 2 (32) 

And 

£ = 0.35k ^ (33) 

dy 

RESULTS AND DISCUSSION 

Calculations are presented for turbulence levels of 1.4 and 2.8 percent at zero pressure gradient and 2.2 percent for a 
favorable pressure gradient with an acceleration parameter (k) of 0.75xl0" 6 . It was necessary, to achieve good 
numerical results, to use about 100 grid points in the cross-stream directions. Comparisons of the calculations were 
made with the experimental data of Blair and Werle (1980, 1981). These experiments were conducted in a low speed 
wind tunnel (30 ft/sec to 130 ft/sec) under ambient air conditions utilizing grids for turbulence generation. A 
favorable pressure gradient was produced by the use of wedge in the test section. The test section was heated at a 
rate of 0.078 Btu/ft 2 -sec with an unheated length of 0.141 ft. The experimental boundary conditions used in the 
calculations are listed in Table IV. 

The resulting calculations using Eqs. (1) to (6) are compared to Blair and Werle’s experimental data in Figs. 6 to 8. 
Figure 6(a) is an example of the abrupt onset of transition that was obtained when intermittency was not used. The 
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TABLE IV— BOUNDARY CONDITIONS FOR BLAIR-WERLE DATA 


Parameter 

Grid 1 

Grid 2 

k = 0.75x10 ^ 
Grid 2 

Free-stream velocity, U e , ft/s 

100.0 

100.0 

32.6 (x = 0) 

Free-stream temperature, T e , °F 

68.5 

68.5 

68.5 

Wall heat flux, q w . Btu/ft 2 /s 

0.078 

0.078 

0.078 

Free-stream turbulence intensity, Tu e , percent 

1.4 

2.8 

2.2 

Free-stream turbulence kinetic energy, k e , ft 2 /s 2 

2.94 

11.8 

0.79 

Free-stream dissipation rate, £ e , ft 2 /s 2 

127 

511 

1.51 


value of using intermittency is clearly demonstrated in Fig. 6(b). There is generally good agreement between 
experiment and prediction. Note in Fig. 7 that even after spot formation begins (transition onset) the boundary layer 
acts as if it were a laminar boundary layer up to intermittency values of 0.40 (1 .4-percent case) and 0.49 (2.8-percent 
case). This is consistent with the measured velocity profiles of Sohn, Reshotko, and O’Brien (1989) which showed a 
Blasius profile for the laminar zone and a laminar like profile for the overall profile, for intermittency values up to 
0.34 with a 1 -percent free-stream turbulence. The numerical calculations capture the position of the minimum heat 



(a) Without intermittency, Tu e = 2.8 percent. 



(b) With intermittency, Tu e = 2.8 percent. 
(«= 11°, P = 0.65). 


Figure 6. — Use of intermittency to model transition region. 
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(a) Total turbulence intensity, Tu e , 1 .4 percent. (b) Total turbulence intensity, Tu e , 2.8 percent 


Figure 7. — Comparison of prediction with experiment (zero pressure gradient, 
a — 1 1 p = 0.65) 
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Distance from leading edge, x, ft 

Figure 8. — Comparison of prediction with experiment (favorable 
pressure gradient, k = 0.75x1 0 -6 , a = 5°, p = 0.8). Total turbu- 
lence intensity, Tu e , 2.2 percent. 


transfer. This is also true in the case with pressure gradient (Fig. 8). Figure 8 shows that, for the case of pressure 
gradient, there is general agreement between experiment and prediction, although the computed curve appears 
abrupt in the region of the minimum heat transfer. There is clearly a need to improve the model. 


CONCLUDING REMARKS 

A demonstration has been made of the important role played by incorporating the characteristics of the transition 
region to provide improved modeling for predictive codes. A model has been developed for predicting the 
nondimensional spot formation rate (N) in terms of turbulent spot variables. This model allows for a determination 
of the transition length and intermittency in terms of the turbulent spot velocity, turbulent spot spreading angle, and 
the turbulent spot physical dimensions. The model produced values of N which are in acceptable agreement with the 
values obtained from zero pressure gradient experiments. There is a need to further confirm the present approach by 
investigating the behavior of turbulent spots under conditions of pressure gradient and curvature for a determination 
of their relationship to transition length and intermittency. It is expected that the above information will aid the 
transition modeling efforts for the improved prediction of shear stress and heat transfer on turbine blades. 
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